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ABSTRACT 

The direct detection of extrasolar planets by imaging means is limited by the 
large flnx of light from the host star being scattered into the region of interest 
by a variety of processes, inclnding diffraction. Coronagraphs are devices that 
snppress the nndesirable scattering of light cansed by diffraction. In a corona- 
graph the sensitivity limit for high dynamic range is limited by the propagation 
of errors introdnced by the imperfect optical system to the final image. In this 
paper we develop theory and simnlations to nnderstand how snch errors propa¬ 
gate in a coronagraph. We describe the response of classical and band-limited 
Lyot coronagraphs to small and large errors in the placement of the central star, 
and identify ways of making snch coronagraphs more robnst to small gniding 
errors. We also nncover featnres of the decentered PSF that can lead to spnrions 
detection of companions, especially with aggressive, high dynamic range coron¬ 
agraphs dedicated to companion searches aimed at finding extrasolar terrestrial 
or Jovian planets. 
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1. Introduction 


It is undesirable that the effect of edge diffraction from the entrance aperture of a 
telescope results in the scattering of light into regions of great interest for study of the 
circumstellar environment of stars. The purpose of coronagraphs is to select or modify the 
spatial frequency content of the light, to effect suppression of diffracted light in a desired 
manner. Discussion of the theory of diffraction limited stellar coronagraphs is typically 
limited to the on-axis point spread function (PSF), often with the assumption of perfect 
optics. In recent years there has been an explosion of new concepts for coronagraphs, many 
of which can achieve contrasts of 10“^° appropriate for Terrestrial Planet Finder (TPF) 
applications in the absence of phase errors. A key question is the tolerance of coronagraphs 
to the variety of imperfections that might be encountered in the real world. Some previous 
studies have incorporated models of phase errors (e.g. Malbet (1996); Sivaramakrishnan et al. 
(2001); Green & Shaklan (2003)), although these have not focused on delivering insight into 
how the errors propagate to the hnal image, and how to design a more robust coronagraph 
at a conceptual level. 

Here we focus our initial analysis on the propagation of tip-tilt errors in Lyot coron¬ 
agraphs (Lyot 1939). A coronagraph is an instrument that suppresses light in a specihc 
position in image space, and thus has a spatially variable PSF. The connection between 
tip-tilt of the wavefront (or equivalently decenter of the focal plane stop) and the response 
of the hnal image plane is important as an error source, and leads to fundamental insight 
into, and understanding of, the operation of a Lyot coronagraph. 

A hard-edged (binary) Lyot coronagraph is remarkably tolerant of tip-tilt errors, even 
for very small focal spots. This is curious, particularly given that one of the most scien- 
tihcally successful coronagraphs, the Johns Hopkins Adaptive Optics Coronagraph (AOC) 
(Golimowski et al. 1992) responsible for the discovery of the hrst bona hde brown dwarf 
(Nakajima 1994) incorporated a tip-tilt system. Whilst it was envisioned that this would 
improve the performance of the coronagraph, in fact there is little beneht to suppression 
of diffracted light, as shown below. Also surprising is the counter-intuitive result that a 
graded or apodized focal plane spot is less sensitive to small tip-tilt errors than a hard-edged 
coronagraph, despite the fact that more light passes through the partially transmissive stop 
when the star wanders off axis. 
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2. Second order monochromatic coronagraphic theory 

The phase on the telescope aperture is 0(a;), where x = {xi,X2) is the location in 
the aperture, in units of the wavelength of the light (see Figure 1). The corresponding 
aperture illumination function describing the electric held amplitude and relative phase in 
the pupil is = A{x) = A[x) (1 + i0(a;) — 0(a;)^/2 + ...), whose Fourier transform, 
Eb = a{k)*{6{k) + i^{k) — <F(A;)*<F(/c)/2 + ...), is the electric held in the image plane B. 6 is 
the two-dimensional Dirac delta function, and k = (ki, /C 2 ) is the image plane coordinate, in 
radians and * is the convolution operator. Our convention is to change the case of a function 
to indicate its Fourier transform. We multiply the image held Eb by a mask function M{k) 
to model the image plane stop of the coronagraph. The image held immediately after the 
stop is Ec = M{k) Eb- The electric held in the re-imaged pupil following the image plane 
stop. Ebb is the Fourier transform of Ec- We use the fact that the transform of Eb is just 
the aperture illumination function Ea itself: 

Eb) = m{x)*EA 

= m{x)*[A{x) {1 + i(l){x) — 

<P{xf/2 + -)] ( 1 ) 

If the Lyot pupil stop transmission is n{x), the electric held after the Lyot stop is Eb = 
n{x)EB). The transform of this expression is the hnal coronagraphic image held strength 
when the wavefront phasor is expanded as a power series in the phase 0: 

Eb = N{k)*[M{k)EB] 

= N{k)*[M{k) {a{k)*{6{k) + i^{k) - 

^k)*^k)/2 +...))]. ( 2 ) 

Understanding high dynamic range Lyot coronagraphy hinges on understanding the structure 
of the held strength Eb> in the Lyot plane located at D. 


3. Guiding error in a Lyot coronagraph 

The ehect of small tip-tilt errors on a Lyot coronagraph operating on a high Strehl ratio 
image is described by a truncated version of equation (1). The mask function in a Lyot 
coronagraph is best expressed as M{k) = 1 — IF(A;), where W{k) is the ‘image stop shape’ 
function. For a hard-edged stop W{k) = Yi{Dk/s), where s is the image stop diameter in 
units of the resolution of the optical system. If the image plane stop is opaque at its center, 
IF(0) = 1 (which constrains w(x) to have unit area). The FT of the stop function M(k) is 
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m{x) = d(x) — w(x), so the Lyot pupil electric field of a Lyot coronagraph can be expressed 
as 


Fd = [(5(a;) — t(;(a;)]*[74(a;) (1 + i0(a;)— 

H^r/2)] (3) 

for sufficiently small phase errors (he., |0| -C 1) in the pupil. Pure tip-tilt error is described 
by a phase function </)(x) = a ■ x = aiXi + 0 : 2 X 2 (o is in radians per wavelength in pupil 
space). We require that the image displacement be much less than a diffraction width, so 
\a\D 1. Following the method developed in Sivaramakrishnan et al. (2002), and truncating 
our expansion at the second order, we derive an analytical expression for the Lyot pupil field 
(which is typically valid for Strehl ratios of the order of 95% and above (Perrin et al. 2003)): 

Ed = [5{x)—w{x)]*[A{x){l+ia-x — 

{a-xf/2+ ...)]. (4) 

Ed is therefore the sum of a zero order term 

Elo = A{x) - w{x)*A{x), (5) 

a first order term 

Eli = i[oi ■ xA{x) — w{x)*{a ■ xA{x)], (6) 

and a second order term 

El 2 = —\[{a-x)‘^A{x)—w{x)*{{a-x)'^A{x))]. (7) 

The behavior of these three terms is most easily understood by following this analysis in 

the case of a band-limited Lyot coronagraph (Kuchner & Traub 2002). We use a coronagraph 
with an image plane stop shape function which possesses a FT of w{x) = n(xi/e, X 2 /e)/e^, 
where e = D/s (s is of the order of a few to 10, and corresponds to the ‘size’ of the image 
plane stop in units of X/D). This simplifies the analytical calculations and brings out the 
salient features of the manner in which tilt errors propagate through a Lyot coronagraph. 
For a hard-edged focal stop, w{x) is a sine function (see Figures 2 and 3). Once we are 
armed with a theoretical understanding of the expressions in equations (5), (6), and (7), we 
can investigate the response of more common Lyot coronagraph designs to guiding errors 
numerically, and also start to address how pupil apodization affects the way guiding errors 
degrade dynamic range. 

The zero order term is well-understood for Lyot coronagraphs {e.g., Sivaramakrishnan 
et al. (2001) and references therein), and is outlined in Figure 2. 
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3.1. First order tip-tilt leak 

The first order term allows light through only at the edges of an unapodized pupil. 
Such behavior is similar to the zero order term. The leaked light can be suppressed by the 
usual undersizing of the Lyot stop. In order to see why this is true, one must consider the 
value of the convolution of the ‘small’ two-dimensional unit-area top-hat function w{x) with 
the function xA{x), as shown in Figure 3. Let us consider an x-tilt (by setting 0:2 = 0). If 
A{x) = 1 inside the pupil, xA{x) is a flat plane with slope ai passing throngh the origin, and 
containing the X 2 axis. The value of the convolution integral when the top-hat fnnction lies 
entirely within the snpport of the apertnre is simply the Xi valne of the offset. Therefore in 
the interior of the pnpil xA{x) = w{x)*xA{x) = Xi. The electric held fnrther than e = D/s 
from the pupil boundary is zero. 

A hard-edged focal stop resnlts in leakage of light into the interior of the pnpil from the 
wings of the sine fnnction (see Figure 3 ). A graded focal stop has a more compact Fourier 
transform than a hard-edged stop. In this case the interior of the Lyot pnpil, where the held 
remains zero, is larger than that of a hard-edged coronagraph’s Lyot pnpil. This resnlts in 
less sensitivity to tilt error for the same Lyot plane stop geometry at high Strehl ratios, even 
when the tilt errors are large enongh to move the star into regions of the focal stop with 
signiheant transmission. 


3.2. Second order tip-tilt leak through 

We apply similar logic to the second order term. In the special case of a clear pnpil, and 
the same band-limited coronagraph, the Lyot pnpil electric held depends on the diherence 
between x^A^x) and w{x)*{xlA{x)). The convolution integral is no longer the identity 
operator even when the top-hat function lies entirely within the pnpil snpport. There is 
a uniform residnal held strength approximately eqnal to everywhere in the interior. 

There is also the same “bright edge” ehect as is seen in the zero and hrst order terms, but 
that is removed by the optimally undersized Lyot stop. The nniform backgronnd in the 
pnpil plane from the second order contribution of a pure tilt term causes a “ghostly PSF” 
to form on axis (not displaced) even with an optimized Lyot stop (see e.g. Figure 5). The 
energy in this PSF varies as the fourth power of the (small) tilt error, and inversely as the 
fourth power of the focal plane stop diameter. First order ehects of defoens will ahect the 
coronagraph in a similar way. It is the combination of these “ghostly PSFs” with the real 
PSF of the star that results in the distorted images shown in Fignre 5. 
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4. The Point-Spread Functions of a Lyot Coronagraph 

Up to this point we have concerned ourselves with small (<C D) tip-tilt errors in Lyot 
coronagraphs. Here we lift that constraint, and examine the morphology of the PSF of a 
Lyot coronagraph over a wide range of stellar locations relative to the spot center. 

We simulated the PSF of a coronagraph when a star is offset from the center of the stop. 
These PSFs are illustrated in Figure 5 using a spot 8A/D in diameter, although we studied 
both smaller and larger stops. We found markedly different morphologies in three regimes. 
When the star behind the spot is displaced a small amount, the PSF looks similar to that of 
the perfectly aligned coronagraph. The rows in Figure 5 show a sequence of locations of the 
central star, beginning at the very center of the occulting spot, with a Lyot stop diameter 
75% of the entrance aperture diameter. When the star is within \/D of the spot edge, the 
PSF develops outcrops that are not at the location of the star. When the star is located at 
the very edge of the spot, or outside it, the PSF takes on a typical direct image PSF shape. 

The three rows of images in Figure 5 are the PSF in the first focal plane, the Lyot pupil 
plane intensity, and the final coronagraphic PSF and shown in radial prohle. We note the 
appearance of the fake source located about 2\/D from the star in the coronagraphic PSF 
at a misalignment of 2\/D. The manner in which placement errors interact with higher 
order errors, such as spherical aberration, has not been studied yet. This suggests that 
PSF modeling of coronagraphic data should be performed with care to avoid misinterpreting 
structure close to the spot edge in the image (e.g. Krist et ah (1998); Krist (2004)). 

This exercise is relevant to coronagraphy on very high Strehl ratio images, although it 
also has immediate applicability to coronagraphic science carried out today, with the HST 
ACS, for instance, if bright structures were present behind the focal stop but near its edge. 

Figure 4 shows coronagraphic rejection efficiency as a function of tilt error for several 
focal plane stops. Typical coronagraphic reductions of the best current space-based data 
demonstrate that imperfect calibration data and temporal variations in the PSF set the 
limits on dynamic range (Krist et ah 1998), so we avoid using simplistic estimates of dynamic 
range using monochromatic simulations to evaluate the actual effects of tilt errors. We use 
the fraction of transmitted central source light as a metric of coronagraphic performance. 
We dehne the transmittance of a coronagraph to be the integrated light in the hnal focal 
plane, excluding the region inside the focal stop (weighted by the focal stop transmission 
for the Gaussian case). This quantity is directly related to the photon-limited noise, albeit 
qualitatively. We choose this quantity as a metric for the purposes of this paper in preference 
to contrast, as it is independent of the choice of an inner working angle for the coronagraph. 
The transmittance of light is calculated for s = 3, 6, and 9 X/D focal stops, with a hard- 
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edged Lyot stop undersized by 0.5 x 1/s of the pupil radius. A Gaussian apodized focal 
stop is compared to the hard-edged focal stop. Since the Fourier transform of an apodized 
stop is more compact, less light bleeds into the center of the pupil (see Figures 2 and 3). 
It is remarkable that the Gaussian apodized stop is more efficient even in the presence of 
qnite large tilt errors. For example, a 3X/D FWHM Gaussian focal stop snppresses more 
light than a 3\/D diameter hard-edged stop even at 0.7\/D tilt (see Figure 4), despite 
the transmission of the stop being 10% at this radius. The remarkable robustness of the 
classical Lyot coronagraph is apparent in Figures 5 and 7, contrary to the expectation that 
led Golimowski et al. (1992) and Lloyd et al. (2001) to incorporated tip/tilt control systems 
into Lyot coronagraphs. The leakage of light from the central star remains concentrated close 
to the edge of the image of the focal stop nntil the central star gets to within a resolntion 
element of the stop edge. This fact, combined with its ease of manufacture and its broad-band 
performance, makes the Lyot coronagraph interesting even in the era of novel coronagraphic 
designs, which must all be well-understood in terms of tolerance to the variety of errors that 
might exist in real telescopes. 

The comparison of Gaussian and hard-edge coronagraphs on an equal footing is com¬ 
plicated by the dehnition of an appropriate eqnivalent width for the Gaussian stop, and the 
nndersizing of the Lyot stop. For the purposes of comparison, we characterized the width 
of the Gaussian stop by a where the transmission of the stop is 1 — exp{—x‘^/2a‘^). We 
adopt the convention of Sivaramakrishnan et al. (2001) and dehne a Lyot tnning paramater 
T which dehnes the fractional radial undersizing of the Lyot stop in nnits oiDjs (or D/a). 
For a hard edges Lyot coronagraph, 7F ~ 0.5 results in most of the performance benehts 
of undersizing the Lyot stop, as the Lyot stop excludes the core of the w{x) sine function 
around the edge of the pupil. Further undersizing in this case results in relatively small gains 
as the wings of a sine function decay slowly (this is calculated in detail in Makidon et al. 
(2000)). For a Gaussian stop, however, the wings are suppressed, and gains continne with 
further undersizing (see Fignre 8). The ultimate application of snch tapering of the focal 
stop to achieve the most compact w{x) is the generalization to more arbitrary fnnetions with 
the concept of the Band Limited Goronagraph (Kuchner & Traub 2002). The rejection of 
snch coronagraphs continues to improve with extremely aggressive undersizing of the Lyot 
stop (see Figure 8. To achieve the very high contrast reqnired for terrestrial planet detection 
a band-limited or gaussian coronagraph with an aggressive Lyot stop (JF > 1.5) For snch a 
coronagraph (see Fignre 6), the near complete rejection of on-axis light is lost with even a 
small tilt error, bnt the coronagraph remains robust against tip/tilt errors in the sense that 
the wings of the PSF are suppressed even for tip/tilt errors of a few \/D. 
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Fig. 1.— The essential planes and stops in a coronagraph. The entrance aperture is A, the 
direct image at B falls on a mask whose transmission function is M{k). The re-imaged pupil 
plane D, after being modihed by passage through a Lyot stop with a transmission function 
n{x), is sent to the coronagraphic image at F. 




Fig. 2.— One dimensional representation of a perfectly aligned hard-edged Lyot corona- 
graph. A band limited stop with a top-hat function bandpass does not have the ringing in 
the wings of the sine function. There is no fundamental difference between these designs 
for the purposes here, since w{x) has approximately the same spatial scale for both. Com¬ 
pared to a hard-edged stop, apodizing the focal stop reduces the ringing in the sine function, 
resulting in less light bleeding into the center of the pupil. 
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Fig. 3.— One dimensional representation of the first-order leak due to tilt error in a Lyot 
coronagraph (see equation (6). The effect of tilt is largely conhned to the edge of the pupil, 
which is already suppressed by an optimized Lyot stop. 
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Coronagraph Transmittance 



I’ilt (?u'D) 

Fig. 4.— Coronagraph transmittance as a fnnction of tilt for three hard-edged and Gaussian 
focal stop diameters. Transmittance is dehned as the fraction of light entering the system 
aperture that propagates to the hnal image plane outside the image of the focal stop. It is 
the total fraction of light that the coronagraph suppresses, not the on-axis null depth. Note 
that the suppression of the 6 A/D hard-edge coronagraph improves with small tilt errors as 
a result of the phasing of the dark/bright Airy pattern with respect to the stop edge. The 
rejection factor asymptotes to the fractional throughput of the T = 0.5 undersized Lyot stop 
(of diameter (1 — \js)D for an sX/D focal stop), which always blocks a fraction of the light. 
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Fig. 5.— Focal, Lyot plane and final coronagraphic image intensities for a hard-edged Lyot 
coronagraph with varying degrees of tilt. The occnlting spot is 8X/D diameter, and outlined 
in red in the focal plane images. The outline of the pupil and undersized Lyot stop are 
shown in red in the Lyot plane images. The outline of the image of the focal stop is shown 
in red in the final coronagraphic image. In the final coronagraphic image, the position of 
the star in the image plane is marked with a cross. If the star is behind the focal stop, the 
peak in the coronagraphic image does not correspond to the position of the star, leading to 
‘fake sources’. The radial profiles show the range from the mean to maximum intensity in 
an annulus centered on the center of the focal stop. The non-coronagraphic Airy pattern is 
shown for comparison. 
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Fig. 6.— Focal, Lyot plane and final coronagraphic image intensities for a Ganssian Lyot 
coronagrapli with varying degrees of tilt. The occnlting mask is cr = 8X/D diameter, and 
ontlined in red in the focal plane images. The ontline of the pnpil and nndersized Lyot stop 
are shown in red in the Lyot plane images. The Lyot stop is nndersized more aggressively 
than in the case of a hard-edged coronagrapli (see Fignre 8 and discnssion in the text). The 
ontline of the image of the focal stop is shown in red in the final coronagraphic image. In 
the final coronagraphic image, the position of the star in the image plane is marked with 
a cross. The radial profiles show the range from the mean to maximnm intensity in an 
annnlns centered on the center of the focal stop. The non-coronagraphic Airy pattern and 
the transmission profile of the mask are shown for comparison. Note that the snppression of 
diffraction in the wings is superior to the hard-edged case, even for tilt errors up to 3X/D, 
at which the Gaussian stop transmission is 7%. 
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Fig. 7.— Combined effect of Lyot stop tuning parameter JF and tilt errors on a hard- 
edged Lyot corongraph indicated by total coronagraph transmittance. A family of Lyot 
coronagraphs with 8 X/D diameter focal stop and varying Lyot stop diameters is shown. 
The no Lyot stop case accounts only for the fraction of energy suppressed by the focal stop. 
The progressive undersizing of the Lyot stop from JF = 0 (a Lyot stop that is the exact 
image of the input pupil) in steps of JF = 0.25 rejects both on-axis and off-axis light. The 
point of diminishing returns is at JF ~ 0.5 as found by Sivaramakrishnan et ah (2001). The 
transmittance asymptotes to the transmission of the Lyot stop. 
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Fig. 8.— Combined effect of Lyot stop tuning parameter JF and tilt errors on a Gaussian 
Lyot corongraph indicated by total coronagrapli transmittance. A family of Gaussian Lyot 
coronagraphs with a = 8X/D focal stop and varying Lyot stop diameters is shown. The no 
Lyot stop case accounts only for the fraction of energy suppressed by the focal stop. The 
progressive undersizing of the Lyot stop from JF = 0 (a Lyot stop that is the exact image 
of the input pupil) in steps of JF = 0.25 rejects both on-axis and off-axis light. Unlike the 
hard-edged case, the on-axis rejection continues to improve to JF > 1 since Gaussian wings of 
a Gaussian continue to drop rapidly unlike the wings of the sincfunction. The transmittance 
asymptotes to the transmission of the Lyot stop. 




